Dynamics of vortex penetration, jumpwise instabilities and nonlinear surface 
resistance of type-II superconductors in strong rf fields 
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We consider nonlinear dynamics of a single vortex in a superconductor in a strong rf magnetic 
field Bo sinojt. Using the London theory, we calculate the dissipated power Q{Bo,uj), and the tran- 
sient time scales of vortex motion for the linear Bardeen-Stephen viscous drag force, which results 
in unphysically high vortex velocities during vortex penetration through the oscillating surface bar- 
rier. It is shown that penetration of a single vortex through the ac surface barrier always involves 
penetration of an antivortex and the subsequent annihilation of the vortex antivortex pairs. Using 
the nonlinear Larkin-Ovchinnikov (LO) viscous drag force at higher vortex velocities v{t) results in 
a jump-wise vortex penetration through the surface barrier and a significant increase of the dissi- 
pated power. We calculate the effect of dissipation on nonlinear vortex viscosity ri{v) and the rf 
vortex dynamics and show that it can also result in the LO-type behavior, instabilities, and thermal 
localization of penetrating vortex channels. We propose a thermal feedback model of which 
not only results in the LO dependence of ri{v) for a steady-state motion, but also takes into account 
retardation of temperature field around rapidly accelerating vortex, and a long-range interaction 
with the surface. We also address the effect of pinning on the nonlinear rf vortex dynamics and the 
effect of trapped magnetic flux on the surface resistance Rs calculated as a function or rf frequency 
and field. It is shown that trapped flux can result in a temperature-independent residual resistance 
Ri at low T, and a hysteretic low-field dependence of Ri{Bo), which can decrease as -Bo is increased, 
reaching a minimum at Bo much smaller than the thermodynamic critical field Be- We propose 
that cycling of rf field can reduce Ri due to rf annealing of magnetic flux which is pumped out by 
rf field from the thin surface layer of the order of the London penetration depth. 

PACS numbers: 74.25. Nf, 74.25. Op, 74.25. Qt 



I. INTRODUCTION 

The behavior of superconductors in strong rf fields in- 
volves many complex mechanisms related to a nonlinear 
electromagnetic response of nonequilibrium quasiparti- 
cles, pairbreaking suppression of superconducting gap A 
and penetration of vortices at higher rf amplitudesi'^ . 
The physics behind the nonlinear rf response has recently 
attracted much attention due to the development of a 
new generation of high-performance superconducting Nb 
cavities for particle accelerators, in which the peak sur- 
face GHz fields B{t) = Bosinujt close to the thermody- 
namic critical field Be were reached at a very high quality 
factor ~ 10^ — 10^"'^ characteristic of the Meissner state^i^. 
At such strong rf fields the peak surface current density 
Bo/fiQ^ approaches the depairing current density at 
which the Meissner state becomes unstable with respect 
to avalanche penetration of vortices once the instanta- 
neous rf field B{t) — Bq sin cut exceeds the superheating 
field Bs ~ Be- In turn, penetration of vortices causes a 
sharp increase in the surface resistance Rg. 

As far as the very high quality factors are concerned, of 
particular interest is the behavior of Rs in s-wave super- 
conductors at low temperatures T <ti Tc and frequen- 
cies w ^ A, for which the rf field cannot break the 
Cooper pairs, and the very low Meissner surface resis- 
tance Rs oc (w^A/r) exp(— A/r) is due to an exponen- 
tially small density of thermally-activated quasiparticles 
(unlike the power-law dependence Rs{T) ~ Ri + CT" due 



to nodal quasiparticles in d-wave superconductors^i^iii^). 
In this case penetration of even few vortices driven by 
extremely high rf currents densities J Jd can pro- 
duce strong energy dissipation comparable to that in the 
Meissner state, which in turn, can trigger thermomag- 
netic flux avalanches and the superconductivity break- 
down. It is therefore important to understand mecha- 
nisms, which control dynamics of single vortex penetra- 
tion under strong rf fields. Yet, the rf field onset of vortex 
penetration B^, and the dissipated power Q as functions 
of Bq and uj, and the relation between By and the ther- 
modynamic Be and the lower critical field Bd are still 
not well understood. These problems include complex ki- 
netics of the emergence of the vortex core at the surface, 
and the subsequent nonlinear large-amplitude oscillation 
of the vortex at the surface driven by strong rf currents 
much higher than depinning critical current density. This 
situation cannot be described by a well developed linear 
electrodynamics of a superconductor in the pinned mixed 
state weakly deformed by rf current a^'^°i^^'^^ . Some is- 
sues of nonlinear vortex dynamics in ramping magnetic 
fields have been addressed in extensive numerical simu- 
lations of the time-dependent Ginzburg-Landau (TDGL) 
equationsi^iiiii^ii^ valid at T « Tc, and molecular dy- 
namics simulations^^. However, few experimental and 
theoretical results on vortices driven by very strong rf 
currents at low temperatures have been published in the 
literature. 

In this paper we address nonlinear rf dynamics of a sin- 



gle vortex moving in and out of a type-II superconductor 
through an oscillating magnetic surface barrier locally 
weakened by a surface defect. We show that already this 
basic problem involves a rich physics, since even weak 
Meissner fields Bq ^ Be can drive the vortex with ve- 
locities v{t) so high that the linear Bardeen-Stephen vis- 
cous drag model becomes inadequate. As a result, the 
vortex velocity v{t) can exceed the sound velocity, caus- 
ing the Cherenkov generation of hypersoun d^^'^^ . More- 
over, v(t) can exceed the critical velocity Vq, above which 
the vortex drag coefficient r]{v) decreases as v increases, 
and the viscous drag force = vri{v) reaches maxi- 
mum at the critical velocity vq, resulting in the jump- 
wise Larkin-Ovchinnikov (LO) instability2£i2i. The LO 
instability has been extensively investigated by dc trans- 
port measurement o^^i^'^i^'^i^^i^^i^''i^^ on both low-Tc and 
high-Tc superconductors for which — 1 ^ 10 kni/s have 
been typically observed at low T and B. Single- vortex 
dynamics under strong rf field also involves annihilation 
of vortex-antivortex pairs, and a cascade of single, double 
and multiple vortex penetrations. Competition of the rf 
driving force, image attraction to the surface, and the vis- 
cous drag force results in a strong dependence of the dis- 
sipated power Q on the rf amplitude and frequency. Very 
high vortex velocities achieved at fields Bq ~ Be required 
to break the surface barrier make it possible to probe 
the behavior of vortices under extreme conditions, for 
which the Lorcntz driving force approaches its ultimate 
depairing limit. Because of strong heating effects, these 
conditions are hard to reproduce in dc transport exper- 
iments (except in high-power pulse measurements22i^) . 
The importance of heating effects for transport insta- 
bilities in superconductors at low temperatures is well 
known^ii^i^i^^. In this paper we show that heating is 
a key limiting factor for the high-field surface resistance 
at T ^ Tc as well, even for single vortices driven by 
strong rf Meissner currents. In particular, viscous vortex 
dynamics coupled with electron overheating can result 
in the LO-type behavior of rjiv), thermal rf breakdown, 
long-range interaction (on scales much greater than the 
London penetration depth) between a vortex and the sur- 
face and between vortices themselves. 

The paper is organized as follows. In section II we 
establish the main parameters of interest by consider- 
ing penetration and dissipation of a single vortex over 
the oscillating surface barrier in type-II superconductors 
described by the dynamic equation, in which the linear 
Bardeen-Stephen viscous drag force is balanced by the 
Lorentz driving force and the image attraction force at 
the surface in the London theory. Even in this basic 
model rf vortex dynamics always involves annihilation of 
vortex-antivortex pairs for Bq w By close to the penetra- 
tion field By and a strong dependence of the dissipated 
power Q[Bq,ui) on the rf frequency and amplitude. In 
section-Ill we show that the Bardeen-Stephen model ac- 
tually has a very limited applicability because vortices 
breaking through the surface barrier reach supersonic ve- 
locity so the dependence of the viscous drag coefficient 
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FIG. 1: Vortex (shown as the filled circle) penetrating by the 
distance u{t) from the semi-infinite surface {x = 0) exposed 
to a uniform parallel rf field B{t). The open circle shows the 
position of the antivortex image. 



rj(v) must be taken into account. In this case vortex 
dynamics becomes strongly coupled with noncquilibrium 
overheating of the vortex core, resulting in jumpwise pen- 
etration of single vortices through the surface barrier, 
and significant increase in Q. In section IV we consider 
the effect of pinning on rf surface resistance. In partic- 
ular, we show that trapped vortices can result in a tem- 
perature independent, field- hysteretic residual resistance, 
which can decrease as the rf field increases. Pinned vor- 
tices can also produce hotspots, which ignite thermal rf 
breakdown. Section V is devoted to dissipation around 
hotspots and their nonlinear contribution to the global 
surface resistance. The thermal breakdown of the Meiss- 
ner state ignited by vortex hotspots is addressed. Section 
VI concludes with a discussion of the results. 



II. PENETRATION OF A VORTEX OVER THE 
OSCILLATING SURFACE BARRIER 

A. Dynamic equations and time scales 

Penetration of vortices in a superconductor is con- 
trolled by the Bean-Livingston surface barrier, which re- 
sults from a competition between the Meissner screening 
currents pushing the vortex in a superconductor and the 
attraction force between a vortex and the surface'^ This 
surface barrier oscillates under the rf field, so motion of 
a vortex in and out of a superconductor is described by a 
dynamic equation. We consider here a type-II supercon- 
ductor in the London theory, assuming that the rf field 
B[t) = Bq sin Lut of amplitude Bq and frequency tu is ap- 
plied parallel to the flat surface at x > as shown in Fig. 
1. Then the equation of motion for a single vortex driven 
by the rf Meissner balanced by the image attraction force 
and the viscous drag force takes the form 



TJqU 



Mo A 



sin ujt~ 



2T:fiQX^ 



, (1) 



where u(t) is the distance of the vortex core from the sur- 
face, A is the London penetration depth, 770 — 4'oBc2/Pn 
is the Bardeen-Stephen vortex viscosity, p„ is the nor- 
mal state resistivity, (j)Q is the magnetic flux quantum. 
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FIG. 2: Snapshots of expanding vortex semi-loop emerging 
from a surface defect (black dot). The quicker expansion of 
the loop along the surface is due to the gradient of the Meiss- 
ner current J{x) oc exp(— it/A) and the LO instability. 



Bc2 = 0o/27r^^ is the upper critical field, Ki{x) is the 
modified Bessel function. Here we introduce the local co- 
herence length at the point of the vortex entry, which 
provides the cutoff in the London theory. For u < ^g, 
the last term in Eq. ([1]) gives a constant force of vor- 
tex attraction to the surface due to the formation of a 
"core string" of depressed order parameter revealed by 
computer simulations of the GL equationsi^i^. 

In this work we treat the emergence of the vortex phe- 
nomenologically, assuming that it first appears in a small 
defect region at the surface. For the results presented be- 
low, the actual nature of the defect is not important as 
long as the defect size is much smaller than A, and the 
local is larger than the bulk coherence length ^. The 
vortex penetrates at the field B{t) > for which the lo- 
cal surface barrier disappears because the peak Meissner 
force (poBo/ fJ-oX exceeds the maximum attraction force 
to the surface (j)lKi{2^s / X) /2fioX^ ■ For <C A, we can 
expand Ki{x) ~ l/x, and obtain 



(2) 



which basically defines in terms of the observed local 
penetration field iJ^,, which has been calculated for dif- 
ferent types of surface defects^i^i^li^. We assume that 
there is a distribution of sparse small regions with re- 
duced local By on the surface where vortices first enter. 
Penetration of straight vortices can only be initiated by 
linear defects (for example, dislocations or grain bound- 
aries) parallel to the vortex line. For more common 3D 
surface defects, such as precipitates or local variation of 
chemical composition, a vortex first emerges as a semi- 
loop, which then expands as illustrated by Fig. 2. The 
initial penetration of a curved vortex in Fig. 2 can hardly 
be described by Eq. ([T]) for a straight vortex parallel 
to the surface. However, the circular vortex semi-loop 
very quickly evolves into a loop strongly elongated along 
the surface because of the gradient of the Meissner cur- 
rent J{x) = {Bq/ hqX) exp{—u/ X) and the LO instability, 
which effectively straightens the vortex due to jump- wise 
lateral propagation of the loop, as shown below. Thus, 
Eq. Il]) can be used after a short transient time, which 
is still much smaller than the rf period. 

Eq. ^ gives By close to the superheating field Bg — 
0.745i?c, above which the Meissner state in extreme type- 
II superconductors with k = A/^ 3> 1 becomes absolutely 
unstable with respect to weak periodic perturbations of 



the order paramete r'^^i'^^i'^-'^i'^^i^'^i'*^ as the Meissner cur- 
rent density at the surface By/ ^qX exceeds the local de- 
pairing current density J^. For Bq > By, a vortex moves 
in and out the superconductor under the action of the 
rf field. Since Eq. ^ can only be used on the scales 
u{t) > we neglect here a possible dependence of 770 
on w, although this dependence can occur if a long-range 
interaction of the vortex core with the surface due to 
nonequilibrium effects is taken into account, as shown 
below. 

To estimate the scale of vortex oscillations and max- 
imum velocities, we first disregard the image attraction 
force, which becomes negligible at distances m > A. Then 
the solution of Eq. ([T]) takes the form 



u{t) ^ Xln 
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(3) 



where to — sin^^ {By /Bq)/lj is the time of vortex entry. 
The maximum vortex penetration depth u„ corresponds 
to cosut = — 1, whence 
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Here Um depends logarithmically on the rf field and fre- 
quency. From Eq. ([3]) we estimate the time r for the 
vortex to move by the distance ~ A from the surface. 
For GHz frequencies and the materials parameters of Nb 
and NbaSn, t turns out to be much shorter than the rf 
period so cosw(to + t) ~ cos ujIq — wrsinojioi hence. 



4>QBy pnS, 



(5) 



Taking p„ = 0.2^flm, A = 90 nm, ^ = 3 nm for NbgSn^^, 
we obtain t ~ 3 x 10^^^ s and wr ~ 0.04 for 2GHz. 
Likewise, taking p„ = lO^^fim and A = ^ = 40 nm for 
Nb, yields r = 4 x lO^^^s. Eq. ^ gives the lower Hmit 
for r because the image vortex attraction increases r. 

Now let us consider rf vortex dynamics in more detail. 
During the positive rf half-period, the vortex penetration 
starts once B{t) exceeds the local By. Because the vor- 
tex currents flow antiparallel to the Meissner currents at 
the surface, penetration of the vortex suppresses the lo- 
cal pairbreaking instability. At the time when the rf field 
almost changes sign, the vortex penetrates by the maxi- 
mum distance u,„ then it turns around and starts coming 
back. However, for negative B{t), the current density at 
the surface J(0, t) is now a sum of the vortex currents and 
the parallel Meissner rf currents. As a result, when the 
outgoing vortex reaches the critical distance Uc from the 
surface, J(0,t) exceeds J^, causing penetration of an an- 
tivortex before the vortex exits. The antivortex is driven 
into the sample by the Meissner current and by the at- 
traction to the outgoing vortex, to which it annihilates 
at the distance Ua from the surface. After that the neg- 
ative B{t) reaches - By, and a new antivortex penetrates 
the sample in the same way as the vortex did for the 
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FIG. 3: Dynamics of vortex penetration and exit calculated 
from Eq. ([T} for t < and Eqs. and dH]) for t > tc for lot = 
0.325 and Bo = 1.02B„. The solid and dashed curves show 
the trajectories of the vortex and antivortex, respectively. 



positive cycle, except that once the antivortex reaches Uc 
on the way out, it creates a vortex at the surface, both 
annihilating at Ua- This process repeats periodically. 

Eq. ID) therefore describes vortex penetration and exit 
until the Meissner current density plus the current den- 
sity of the outgoing vortex being at x = Uc reaches - 
By/ hqX at the surface at the time t — tc defined by: 



Bo|sinwt,| + -^Xi(^) =By 



(6) 



The second term in the l.h.s. of Eq. ^ is twice the 
current density of the vortex at the distance Uc in an 
infinite sample because the outgoing vortex and its an- 
tivortex image contribute equally to J(0,t). For t > tc, 
the vortex with the coordinate and the antivortex 

with the coordinate U-{t) move toward each other, as 
described by the following equations: 
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FIG. 4: Dependence of the second vortex penetration field 
B2 on the dimensionless frequency n = lot / k. 



Um — 3A, the critical distance is Uc ~ 2A, and the vortex- 
antivortex annihilation occurs at Ua — X. For lot <C 1, 
the vortex first accelerates rapidly, penetrating by the 
distance ~ A during a time ~ r, and then slowly turns 
around during the time of the order of the rf period and 
annihilates in a short time ^ t. 

The above results are limited to the field region B^ < 
Bq < Bs where single vortices penetrate independently 
through regions where the Bean-Livingston barrier is lo- 
cally suppressed by surface defects separated by distances 
> A. The case Bq > Bs corresponds to a global pair- 
breaking instability causing multi-vortex avalanche pen- 
etration. Yet, even for B^ < Bq < Bg, a, multi- vortex 
chain penetration is possible. Indeed, penetration of a 
single vortex for B{t) > By suppresses the local pair- 
breaking instability at x = 0. However, as B{t) increases, 
the Meissner current density increases, while the counter- 
flow of surface current density at a; = from the vortex 
decreases as it moves further away from the surface. As 
a result, J{0,t) can again reach Jd, causing penetration 
of the second vortex at t — t2 when the first vortex is 
located at a; = U2 The condition of the second vortex 
penetration is similar to Eq. ([6]), 



V A 



X 



B2 sin ujt2 



(9) 



(8) 



These equations reflect the balance of interaction forces 
between the vortex and antivortex and their correspond- 
ing images similar to those in Fig. 1. The first term 
in the r.h.s. of Eq. ([8]) has the minus sign because the 
Meissner current drives the antivortex in the opposite di- 
rection as compared to the vortex. The initial conditions 
for Eqs. ([7]) and (O are: u+(tc) = "c, u^{tc) = 0; and 
the condition u+{ta) = U-(ta) = Ua defines the annihila- 
tion distance Ua and time ta- 

Dynamics of vortex penetration and annihilation is il- 
lustrated by Fig. 3 where the vortex penetration depth is 



except for the minus sign in the l.h.s. Eqs. ([T]) and ^ 
define the critical rf amplitude B2 below which only the 
single vortex penetration occurs. Shown in Fig. 4 is the 
curve B2{lo) obtained by the numerical solution of Eqs. 
$1]) and ([9]) for NbaSn. These results can be described 
well by the power law dependence 



B2 = [l+ pn°']By, n = 2fiQX^oj/p„ 



(10) 



where = lot/k, a = 0.73 and p — 0.23. For 57 ^ 1, 
the field B2 is close to By, however, dissipation produced 
by penetrating vortices can significantly reduce both By 
and B2 (see below). 
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The field region of the single- vortex penetration _B„ < 
Bo < B2 defined by Eq. pH)) shrinks as the frequency 
decreases. In this narrow field region the effect of vortex 
viscosity can radically change the dependence of Q on 
Bq and Lu. Shown in Fig. 5 are the results of numerical 
solution of Eqs. ([T])- ([5]) and (fTT|) for different frequen- 
cies and the GL parameters k = X/S,s- These data are 
described well by the following formula 



Q = 



TT^O 



Bt 



Bl 



T2 



which reduces to Eq. (fTH]) for u ^ Q. As follows from Fig. 
5 and Eq. (fT4|) . the power Q decreases as oj increases be- 
cause of retardation effects due to vortex viscosity during 
the short fraction of the rf period in which B{t) > By. 



FIG. 5: Dissipated power Q as a function of the rf amplitude 
in the region of the single-vortex penetration Bv < Bq < B2 
for different frequencies ujt: 0.061 (1); 0.162 (2); 0.325 (3); 
0.65 (4). Here Qo = lujcjmB^ /tx^o. 



B. Vortex dissipation 

The power Q = {uir] / 2tt) f dt dissipated due to the 
work of the viscous drag forces is given by 



Q = 



LOT] 



to 



_)dt 



(11) 



where u{t) is the solution of Eq. (IT|), which describes 
dynamics of a single vortex driven by rf field until t — 
tc when the antivortex appears. The second integral in 
Eq. ([TT|) is due to the collapse of the vortex-antivortex 
pair described by Eqs. ([7]) and ([5]). For a quasistatic 
field, Q can be obtained from the change of the vortex 
thermodynamic potential G{u) 
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Mo 
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where i?ci — eMo/0o is the lower critical field and e is 
the vortex self-energy. If the ac field B{t) varies very 
slowly {n <C 1), the dissipated energy equals the sum of 
AG+ — G(0) — G{um) during the positive half-cycle and 
AG- = G{um) — G{Q) during the negative half-cycle. For 
any closed vortex trajectory, which starts and ends at the 
surface, contributions to Q due to vortex self-energy and 
the work § Fi{u)udt of the potential image force Fi{u) 
vanish. Thus, Q is only determined by the work of the 
driving Lorcntz force, ~ (20o^i;//^o)[l ~ exp(— m,„/A)], 
for both vortex and antivortex cycles, where we took ac- 
count of the fact that the main contribution to Q comes 
from the initial acceleration of the vortex during the 
time ^ T when the field B{t) is close to B^. Neglect- 
ing exp(— Mm/A) ^ 1, we have 



(13) 



III. INSTABILITIES AT HIGH RF FIELDS 

Once the field B{t) exceeds B^, the vortex rapidly ac- 
celerates reaching the maximum velocity Vm — A/r: 



(15) 



Eq. (fT5|l gives Vm — 30 km/s for NbaSn and Vm — 10 
km/s for Nb. Not only arc the so-obtained values of Vm 
much higher than the velocity of sound, they may even 
exceed the critical BCS pairbreaking velocity^, 



A 



(16) 



where ^ = hvp/T^A, is the Fermi velocity and m is 
the electron mass. Indeed, taking ^ = 40nm and the free 
electron mass m, we obtain wa =0.8 km/s < Vm for Nb, 
and va = 10 km/s < Vm for ^ = 3nm in NbaSn. Here we 
use the Bardeen-Stephen model for qualitative estimates 
only, ignoring many still not well understood mechanisms 
essential at low temperatures, for example the effect of 
quantized electron states in the core and the core shrink- 
age due to the Kramer-Pesch effect^, resulting in the 
factor ~ ln(Tc/r) in the Bardeen-Stepnen formula^i^i^. 
Yet, for strong rf fields, Bq ^ Be, the linear viscous drag 
force derived for small vortex velocities, becomes inade- 
quate. It was first predicted theoretically?? and observed 



in many experimenti 
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that the dependence 



of ?] on V at high vortex velocities results in a nonmono- 
tonic viscous drag force /„ = vri(v) and jump-wise insta- 
bilities. 



A. Instabilities of viscous flux flow 

A nonlinear viscous drug force was first calculated 
by Larkin and Ovchinnikov (LO)^, who showed that 
nonequilibrium effects in the vortex core decrease the 
drag coefficient rj as v increases: 



ri{v) 
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(17) 
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FIG. 6: The viscous drag force fv{v) as a function of the 
vortex velocity for the LO instability. The dashed curve shows 
fv{v) at rjo — 8rii. For 770 > 8771, the N-shaped dependence 
fviv) develops, as illustrated by the solid curve plotted for 
rji = 0.05r;o. The arrows show the jumpwise change of v{F), 
as the driving force F increases above the maximum value of 
/„(!) and then decreases below the minimum value /u(3). 



where 770 is the Bardeen-Stephen viscosity. The critical 
velocity vq in the dirty limit is given by: 



vo ~ 0.6 



1/2 



T 



1/4 



(18) 



Here ^i is the mean free path due to impurities, and t^{T) 
is the quasiparticle energy relaxation time. Eq. (fT7|) 
results in a nonmonotonic dependence of viscous drag 
force, = vrj(v) on the vortex velocity: 



riiv 



1 + v'^/v^ 



(19) 



Here, following the LO approach, we use two effective 
viscosities 771 and 772 , where 772 phcnomcnologically takes 
into account the transition to the normal state as v 
reaches the pair breaking velocity ua. For 771 < 8772, 
the force fv(v) always increases as v increases, but for 
rji > 8772 the N-shaped dependence fv{v) develops, as 
shown in Fig. 6. For 772 = 0, and 771 = 770, the drag force 
reaches the maximum value = riQVo/2 at v = vq. 

The LO instability was originally associated with ac- 
celeration of normal quasiparticles in the vortex core by 
electric field, which can increase their energy above A. In 
this case quasiparticles can escape the core if the diffusion 
length Ld = {Dt^)'^/'^ exceeds the core size. The result- 
ing quasiparticle depletion in the core reduces the core 
size and the vortex viscosity according to Eqs. p?|) and 
(fTSll . However, Eq. pT|) is actually more general and may 
result from several different mechanisms. In particular, 
the velocity dependence P?|) can result from coupling of 
the vortex motion with other diffusion process, including 
quasiparticle or temperature diffusion around the mov- 
ing vortex core. For example, the electron overheating of 
the core can lead to Eq. ([T7]) as follows. 



The power 777;^ generated by the viscous drag increases 
the electron temperature in the core T^, reducing the 
vortex viscosity ri{Tm) oc Bc2{Tm)- Linearizing rjiTm) — 
770(1 — Tm/Tc), we write the thermal balance condition 



(T„ - To)g = 7;o(l - T„jTc)v'^ 



(20) 



where g ~ nk/ \n{Lg/ L^) defines the heat flux from the 
core due to the thermal conductivity k(T), where the 
thermal length Lg is of the order of the film thickness 
d for ideal cooling, and ^ is the length related 
to the amplitude of vortex penetration. This logarithmic 
factor will be calculated in the next section in more detail; 
here we just use \n{Lg/L^) ~ In^d/y/^) for qualitative 
estimates. Solving Eq. (|20p for Tm results in the velocity- 
dependent 77(10) of the LO form: 



77 = 770(1 -ro/r,)/(l+77oi;V<?), (21) 
(3/'7o)'/' - [7:kp,JcboB,2iO)HLg/L^)]'^^ (22) 



Substituting here the low-T quasiparticle thermal con- 
ductivity k ~ fc„(A/T)2exp(-A/r), and using the 
Wiedemann-Frantz law fc„p„ — {'!TkB/e)'^T/3, and 
Bc2{0) ^ 0o/27r^o^i and <^o ~ hvp/^ in the dirty limit, 
we reduce Eq. (|22p to Eq. (US)). Here the time constant 



hT (Le\ /A 
_ln(^_jexp ^- 



(23) 



exhibits the exponential temperature dependence similar 
to the energy relaxation time between quasiparticles 
and phonons^ in the LO theory. However, the exponen- 
tial dependence of t^{T) in the thermal model is cut off 
at lower T where k is limited by phonons. 

To estimate vq, we take p„ — 0.2/xr2m, k — 
10"2^/mK, Bc2 = 23T for NbgSn at low temperatures^^ 
and ^ ^ 16 nm, Lg ^ d ^ 1 mm, ln(Le/L^) ~ 
11. For these parameters, Eq. ([22|) gives vq ~ O.lkm/s, 
much smaller than the estimates for wa and Vm- Thus, 
overheating does result in the same Eq. (fT7|) . although 
in this case vortex core expands as it becomes warmer at 
higher velocitiesS^j^, in contrast to the LO core shrink- 
age. Moreover, the critical velocity wo defined by Eq. ([^ 
remains constant at Tc, unlike vanishing wo(2o) for the 
LO mechanism, which dominates at T ~ T^. However, 
for T Tc, both Eq. HID and ^ predict the critical 
velocity vo{T) to exceed the linear viscous drag-limited 
velocity v^iT) cx (1 - T/TcY^^ given by Eq. As a 

result, Eq. ([1]) adequately describes rf vortex dynamics 
at strong fields Bq ~ Be and temperatures close to Tc- 

To evaluate the overheating mechanism in more detail, 
we assume that r]{Tm) depends on a local electron tem- 
perature Tm{t) in the vortex core. The distribution of 
T(r, t) around a moving vortex is described by a thermal 
diffusion equation, 

Cf ^ kV'T~a{T~To)+viT^y{t)f[x-uit),y], (24) 

which, after re-definition of the coefficients, can be re- 
duced to the same mathematical form as the diffusion 
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equation for nonequilibrium quasiparticles. Here C is 
the heat capacity, k is the thermal conductivity, u(t) is 
the coordinate of the vortex core moving with the veloc- 
ity and the function f{x,y) accounts for the finite 
core size, so that j f{r)(Pr = 1. The term a{T — Tq) 
describes heat exchange with the environment. For ex- 
ample, a = h/d in a thin film of thickness d where h is 
the Kapitza conductance at the sample surface, and Tq 
is the bath temperature. For electron overheating, the 
parameter a = C/t^ describes heat exchange between 
electrons and the lattice, where C is the electron specific 
heat and is the time of inelastic scattering of quasi- 
particles in the vortex core on phonons^i^i^i^. The 
last term in the r.h.s. of Eq. ([24|l describes dissipation 
in the vortex core proportional to the viscosity r]{v,Tm) 
taken at the local core temperature T™ (t) , which, in turn, 
depends on v{t). The core form factor /(r) is modeled 
by the Gaussian function /(r) = 7r~-^^f ^ exp(— r^/^^), 
where the core radius can be smaller than ^ at T <C Tj, 
due to the Kramer-Pesch effect^ ( the solutions of Eq. 
(IM)) depend weakly on ^i). We also consider weak over- 
heating for which dependencies of C, k and a on T can 
be neglected. As shown below, T can be regarded as ei- 
ther the electron or the lattice temperature, depending 
on the time scale of the vortex dynamics involved. 

The solution of Eq. given in Appendix A, results 
in the following integral equation for the temperature 
T„i{t) in the vortex core moving with a time-dependent 
velocity v{t) near the surface: 



0.4 r 



Tm{t) = To + 



j_ dt'qjt - ty-^'/'o 

^ Jo 4t' + 1, 

uit) ~ u{t - t')]^ 



{exp 



exp 



[At' + ts)D 
[u{t) + u{t - t')f 
(4f + ts)D 



± 



(25) 



where q(t) — r/[Tm(i), ■i;(i)]u^(i) is the time-dependent 
power generated by the moving vortex, D — k/C is the 
thermal diffusivity, ts = ^1/ D \s the diffusion time across 
the vortex core, and tg = C/a is the electron energy 
relaxation time. The second term in the parenthesis de- 
scribes the effect of the surface: the plus sign corresponds 
to the thermally-insulated surface, dxT{x,t)\x=o — 0, 
and the minus sign corresponds to the ideal cooling, 
T(0,t) = Tq. Here we do not consider microscopic ther- 
mal gradients inside the vortex core^^, assuming that 
their effect is included in the bare 77. 

The integral Eq. (pS)) takes into account retardation 
effects due to diffusive redistribution of r(r, t) around 
an accelerating vortex, so Tm{t) depends on the vortex 
velocity v{t — t') at earlier times. The effect of the surface 
makes T„i(t) dependent on the vortex coordinate u{t) as 
well. Eq. ((25)) simplifies considerably ifv{t) varies slowly 
over the relaxation time tg, and tg, and u(t) > ^1. 
Then q{t ~ t') can be taken out of the integral at t' = 0, 
and Eq. (P5|) yields the following equation for the local 




FIG. 7: Temperature spike in the vortex core after the jump 
described by Eq. (gT)) for tg = lOtd, t^i = Su^ /AD, STo = 
WD/nkSu'^, and Ta = To + riovi\n{Le/i)/2nk. 



temperature difference ST, 



T„ 



To: 



ST,n = 



27rfc 



(26) 



Here I = ^16^/72 « 0.67^i, and 7 = 0.577 is the Euler 
constant^. The second term in the parenthesis decreases 
exponentially for u > Lg, and logarithmically, Ko{z) ~ 
ln(2/z) — 7 for ^ <C u ^ ig. In this case the expression 
in the parenthesis reduces to ln{Lg/£^) +ln{Lg/u), where 
u = ue^ 12. Taking the characteristic amplitude of vortex 
penetration u ^ A, we can present the logarithmic part 
in the form 2\n{Lg/ L^), where = {"^SyY^"^ ^ was 
used before to obtain Eq. (|17p in the thermal model. The 
weak logarithmic dependence of vo on and Lg makes 
Eq. (22) nearly insensitive to the details of heat transfer 
and the behavior of u{t). 

In the other limiting case of very rapid variation of v[t) , 
the vortex reaches the critical velocity vo and then jumps 
by the distance 5u, dissipating the energy W during the 
short time 5t. Then q{t) = qo + WS{t) and Eq. ^ 
results in the following implicit equation for 5Tm{t) at 
t > St: 



xrr u\ '?o'"o , Lg 



WXTrnl^ 

Ankt 



'Su^ /4Dt-t/tg 



(27) 



which describes a temperature spike in the core followed 
by relaxation of 5T„i{t), as shown in Fig. 7. Here the 
first term in the r.h.s. gives ST„i before the jump, and 
the effect of the surface is neglected. 

Next we consider the steady-state temperature field 
T(r) averaged over rf oscillations, where T(r) is deter- 
mined by the balance of the vortex heat source and ther- 
mal diffusion. Solution of the thermal diffusion equa- 
tion in Appendix A yields the following distribution of 
(5r(r) = T(r) — Tq from a heat source localized at the 
thermally-insulated surface (a; = 0) of a slab of thickness 
d, ideally cooled from the other side, ST{d) = 0: 



, cosh^ 

q[x')\n 2iL 

cosh ^ 



' 2d 



'-dx' (28) 
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FIG. 8: Temperature variation (5r(x, J/) — T{x,y)~To around 
a lieat source given by Eqs. pSI) and p9|) for Um ~ 10~'^d/-K 
and normalized by Tq — Q/2nk. 



Here q(x) is the power density averaged over the rf pe- 
riod. On scales greater than the size of the heat source, 
ST{r) depends only on the total power Q — q{x)dx: 



ST(r) - ^ cosh(7ry/2d) + cos(7rx/2d) 
27rfc cosh(7ry/2(i) — co8{nx/2d) ' 



(29) 



Here 5T{r) decays exponentially over the length 2d/7r as 
shown in Fig. 8. Near the heat source 6T{r) weakly 
(logarithmically) depends on details of q{x). As shown 
in Appendix A, the distribution of 6T{0, y) along the 
surface to the logarithmic accuracy is given by 



^my) = -^inV?^ 



« 



(30) 



where ro quantifies a size of dissipation source. The ac- 
count of finite tq cuts off the logarithmic divergence in 
Eq. (j29p at a; = y = 0, resulting in a maximum temper- 
ature disturbance at y = 0: 



o7m ^ — In 



(31) 



Eq. (PT|) reduces to Eq. with u ^ and Lg ^ d. 

The physical meaning of T in the above formulas de- 
pends on the relevant vortex time scales. For example, 
for the supersound vortex penetration or vortex jumps 
on the time scale much shorter than the electron-phonon 
energy relaxation time, the quasiparticle are not in equi- 
librium with the lattice, and T(r, t) in Eq. (P7)) can 
be regarded as an electron core temperature. However, 
steady-state vortex oscillations in the rf field generate 



a dc power, which must be transferred to the coolant 
through phonons. In this case Eqs. ([^ describes 
the lattice temperature distribution around a vortex if 
the phonon mean free path is shorter than the film 
thicknesg^^. Thus, the vortex oscillates in a "warm tun- 
nel" with the lattice temperature 5T{y) shown in Fig. 8, 
but in addition to that the vortex core gets overheated 
with respect to the lattice during short periods of rapid 
acceleration, jumps or annihilation with antivortices, as 
described before. 



B. Jumpwise vortex penetration 

For the LO vortex drag coefhcient 77(17) given by Eq. 
p?)) . the equation of motion becomes 



770U 



1 -I- v?/vl 



-e 



/A 



27r/ioA-^ 



Mo A 



sin uit — 



(32) 



The nonmonotonic velocity dependence of the viscous 
drag force in the l.h.s. of Eq. (|32p qualitatively changes 
vortex dynamics as v(t) exceeds the critical value vq for 
which the viscous force reaches the maximum Fm — 
riQVo/2. Indeed, the differential equation for u{t) has the 
form r]ou/{l + u^/vq) = F{u,t), where F is the net elec- 
tromagnetic force given by the r.h.s. of Eq. ([5^ . We can 
introduce the ratio P of the maximum Lorentz driving 
force at Bq = By to the maximum viscous force: 



P = 



(33) 



As shown above, the Bardeen-Stephen viscous flow re- 
sults in unphysically high vortex velocities, indicating 
that P 3> 1 and the dependence of 77 on u must be 
taken into account. However, in this case there are re- 
gions at the surface where F{u,t) exceeds Fm as shown 
in Fig. 9. In these regions the force balance Eq. ([5^ 
cannot be satisfied and the vortex jumps to the place 
where F{x) < Fm- To see how it happens, we present 
the quadratic equation Eq. ([5^ for ii in the form 



v±{F) 



F{u,t) 



i± Wi- 



F^{u,t) 

~F2 



(34) 



where Fm = r]oVo/2. For v <^ vq and F <C Fm, Eq. 
with the minus sign in the brackets reduces to Eq. ([T]). 

Penetration of the vortex at B{t) > By is therefore de- 
scribed by the first order differential equation u = v- (F) , 
which is well defined only if F{u,t) < Fm, otherwise the 
driving force exceeds the maximum friction force, and 
the square root in Eq. ([34)) becomes imaginary. Vor- 
tex dynamics in this case can be understood from Fig. 
9, which shows an instantaneous profile of F{u, t) for 
B{t) > By. Here a vortex enters the sample with zero 
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u/X 

FIG. 9: Instantaneous profile of the force F{u) normalized 
by Fo = (fioB^/fioX for B{t) = and k = 20. The dashed 

line shows the maximum viscous force F,n = riovo/2, which 
can balance F{u) only in the regions where F{u) < Fm- Thus, 
vortex first moves from u = to u = iti, then jumps from 
point 1 to point 2 after which it moves continuously as de- 
scribed by Eq. p2p until the next jump and annihilation with 
the antivortex on the way back. 
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FIG. 10: Jumpwise vortex penetration calculated by solving 
Eq. ((32]) for k = 34, = 0.244, Bo = 1.02B„ and: «o = 
0.0527;™, P = 38.4 (1), and vq = 0.527J™, P = 3.84 (2). The 
dashed lines show the antivortex penetration. 



velocity at t = to and then accelerates because the net 
force F{u^ t) increases as the vortex moves away from the 
surface and the vortex-image attraction weakens. This 
part of u{t) is described by the equation ii = v-{F) un- 
til the vortex reaches the point 1 where F{u,t) = F,n- 
In the region u > ui, the friction force cannot balance 
the driving force, so the vortex jumps to a point 2 where 
F{u, t) = Fm and the viscous drag can balance the driv- 
ing force. As follows from Eq. (P7|) . the core temperature 
does not change right after the jump {t ^ 0), but then 
starts increasing. After that the smooth parts of u{t) are 
described by to the same dynamic equation ii = v^{F) 
for the vortex, which reaches the maximum penetration 
depth Um then turns around and accelerates toward the 
surface during the negative rf half period. However, as 
the vortex reaches the negative critical velocity —vq on 
the way out, it jumps again and either exits the sample 
or collides with the incoming antivortex and annihilates 
as described in the previous section. The positions of the 
jumps Uj at the corresponding times t = tj are deter- 
mined by the equation F{uj, tj) — ii*"™: 

= ; smwi,-- — ^^Ki[ , (35) 

2 fXoX 27r^oA3 V ^ / 

where the ± signs should be taken for the incoming and 
outgoing parts of respectively. Shown in Fig. 10, 
are examples of vortex penetration dynamics calculated 
numerically from Eq. ()32p . In this case the antivortex 
jumps in and annihilates with the outgoing vortex before 
the vortex reaches the critical velocity —vq. 

As mentioned above, the LO instability facilitates a 
quick evolution of the vortex semi-loop originating at a 
3D surface defect into a straight vortex parallel to the 
surface. Indeed, as evident from Fig. 2, the vortex 
propagation velocity is maximum for the segments of the 
semi-loop perpendicular to the surface because they are 
driven by the continuously increasing maximum Lorentz 



force Fl = B(t)(j)o/ iiqX, while experiencing no counter- 
balancing image forces. As a result, the LO instability 
first occurs for the perpendicular vortex segments, caus- 
ing them to jump along the surface to the place where the 
viscous force is able to balance the Lorentz force. How- 
ever, unlike the parallel vortex segments whose jump dis- 
tance X perpendicular to the surface is limited by the 
London screening, the jump length along the surface of a 
semi-infinite sample is infinite because F^ remains con- 
stant. Thus, the vortex semi-loop turns into a straight 
vortex in a jumpwise manner when the lateral velocity of 
the perpendicular vortex segments reach ±vq. 

Several points should be made regarding the jumpwise 
vortex dynamics. First, the vortex trajectory u(t) com- 
prised of the jumps connected by smooth parts described 
by the equation ii ^ (F) occurs only if F{uj, tj) < at 
the jump points where v{tj) — ±vo- However, at higher 
frequencies, or as the vortex overheating is taken into ac- 
count, there are situations when F{uj,tj) > 0. In this 
case the vortex velocity after the jump exceeds vq and 
the smooth parts of u{t) are described by both branches 
it = v±(F), as shown in Appendix B. 

The second point is that, for the overdamped dynamics 
described by Eq. (j32p . the jumps occur instantaneously 
unless the second ascending branch due to 772 term in Eq. 
(fT9)) is taken into account. However, this branch in the 
LO model corresponds to very high velocities ^ ua, for 
which the adequate theory of the nonequilibrium vortex 
core structure and the vortex drag force is lacking. In our 
phenomenological London approach we assumed that the 
vortex jumps to the nearest point U2 where the friction 
force is able to balance the driving force F{u2) — Fm- 
However, the instantaneous LO dependence 77(1;) does 
not include retardation effects due to finite relaxation 
times of the superconducting order parameter, or dif- 
fusive redistribution of nonequilibrium quasiparticles or 
temperature around a rapidly accelerating/decelerating 
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vortex core. These effects are taken into account by the 
integral Eg. ips)) for the core temperature Tm{t), which 
shows that the vortex jump time and length are affected 
by intrinsic dynamics of 77. Thus, there is a diffusion time 
scale 5t ~ 5v? /D for the vortex jump by the distance 5u, 
where D equals either fc/C in the thermal model or the 
quasiparticle diffusivity in the LO theory. In the thermal 
model this estimate gives 5t ^ A x 10"^^ s if we take 
(5m ~ A = 65nm, k ~ O.lW/mK and C ~ IQOJ/m^K for 
NbsSn at 2K^, or even much shorter time for Nb, for 
which A = 40nm and k ~ lOW/mK. The so-estimated 5t 
is smaller than the inverse gap frequency, indicating that 
once the overheated core gets in the region where the 
friction force is able to balance the driving force, it cools 
down very quickly due to the electron component of ther- 
mal conductivity. In Nb the electron thermal conductiv- 
ity k (X exp(— A/r) is still significant down to T > IK, 
but at lower temperatures 5t may be limited by much 
slower phonon irradiation from the overheated core. At 
the same time, the electron temperature relaxation time 
Te outside the core results from a slow phonon-mediated 
recombination of quasiparticles^ 

T£ ^ ro (^^^ exp , (36) 

which yields 30 ns much longer than the thermal dif- 
fusion time St for Nb at 2K. For lut^ <C 1, the quasiparti- 
cles are overheated with respect to the lattice, in a highly 
inhomogeneous way according to the distribution of the 
lattice temperature T(r) shown in Fig. 8. In this case 
the condition ujT^{r) > 1 of the electron overheating can 
locally be satisfied in colder regions away from the core, 
but near the vortex core the electron temperature can 
be close to the lattice temperature if uJT^{r) < 1 because 
of higher r(r), which greatly accelerates the energy ex- 
change between electron and phonons. For example, for 
Nb at To = 2K, the local increase of the lattice tempera- 
ture to T = 4X yields t,{T) ~ t,{To) exp(-A/To+A/T), 
giving Tg ~ 0.3 ns, and lot^ 1 at 1 — 2 GHz. 

Another contribution to 6t comes from a finite vor- 
tex mass M. In the Suhl model M ~ 2mkp/TT^ is due 
to localized electrons in the vortex core, where kp = 
(Sir'^ny^^ is the Fermi wave vector^'*. The jump by A 
due to the driving force F — (j)QBy/ iiqX takes the time St 
set by the Newton law 2\/St^ ~ F. Using = (/)o/47rAC, 
A^ = ^Qm/ne^, ^ = Hvf/tt^, we obtain 

(37) 

A more accurate account of quantized levels in the vor- 
tex core^^i^ or lattice deformation around the moving 
vortexi^i^'^ can increase the vortex mass, thus further 
increasing the jump time St. Yet, although vortex jumps 
are quantified by the multiple relaxation times discussed 
above, they seem to occur much faster than the rf periods 
we are dealing with in this work. 
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FIG. 11: The dependence of the second penetration field B2 
on the dimensionless frequency — ljt/k, for the LO insta- 
bility at K = 34 and vo = 0.52«m. 

C. Rf dissipation 

Power Q dissipated with the account of the jumpwise 
vortex instabilities can be written in the form 

Q = -{^ i]{u)u^dt + <j> ■q{u){u\ + ul)dt + 

)]}, (38) 

m 

where the integrals are taken over all smooth parts of the 
vortex trajectory u{t), including the vortex-antivortex 
annihilation parts, like in Eq. (jlip . The last term in 
Eq. ([55]) is the sum of energies released during all vortex 
jumps aX t — t„j, from u = u„j_ to u = u„i+, where the 
instantaneous free energy G{u,t) is given by Eq. (fT2|) . 
Shown in Fig. 11 is the second vortex penetration field 
B2 calculated by solving Eqs. and ^ nu- 

merically for the parameters of NbaSn. Here B2 cannot 
be fit with a power-law similar to Eq. (fTO| . 

The LO instability makes the behavior of Q{Bo,uj) 
more complicated as compared to Q{Bq) described by 
Eq. (fT4|) . As shown in Fig. 12, there are three distinct 
field regions of very different vortex dynamics. The pure 
Bardeen-Stephen dynamics like that shown in Fig. 3 is 
limited to a very narrow region of Bq close to By (labeled 
by a in the inset). In this case the vortex penetration 
depth Um turns out to be smaller than ^, indicating that 
the London theory combined with the Bardeen-Stephen 
drag cannot give a self-consistent description of vortex 
dynamics at low temperatures. However, Eq. ([T]) ade- 
quately describes rf vortex dynamics at higher T close 
to Tc where the LO instability is irrelevant because the 
critical velocity uq oc (1 — T/TcY^^ becomes larger than 
v^oc{l- T/T,f/\ 

The parts of the Q{Bq) curves labeled h in Fig. 12 cor- 
respond to an intermediate case, for which the most part 
of u{t), including the initial acceleration of the vortex, 
reaching the maximum penetration depth m„i, and turn- 
ing back, does not involve the LO instability. However, 
as the velocity of the exiting vortex exceeds —vq, it jumps 
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FIG. 12: Q{Bo) calculated with the account of the LO in- 
stability for K = 34, vo — 0.52vm, and 57 = 0.019 (1) and 
n = 0.089 (2). Solid lines show the parts of Q{Bo), for which 
Um > (,s- Inset shows the detailed behavior of Q{Bo) close to 
Bv , where the parts labeled a correspond to the pure Bardeen- 
Stephen dynamics for the entire u{t). Parts b are partially 
controlled by the Bardeen-Stephen drag for the vortex reach- 
ing Um, then turning back and then jumping to the surface 
from u — Uc < Um due to the LO instability. Parts c corre- 
spond to the LO dynamics similar to that in Fig. 10. 



to the surface and disappears, significantly increasing the 
dissipated power Q. Further increase of Bq corresponds 
to the parts of Q{Bq) curves labeled c, for which the LO 
instabilities occur both on the penetration and the exit 
parts of u{t), hke those in Fig. 10. In this case Q{Bq) 
jumps up to a much higher level Q ^ Qq until the sec- 
ond penetration field B2 is reached. Here the behavior of 
Q{Bo) can also depend on uj: for lower frequency (curve 
1), Q{Bo) increases weakly between 0.4 < e{Bo) < 0.8, 
but for e= {Bo- By)/{B2 - S„) > 0.8, the power Q{Bo) 
jumps down. This behavior reflects the change in the 
vortex dynamics: for 0.4 < e < 0.8, the vortex jumps 
out of the sample before the antivortex enters, while for 
e > 0.8, the jump of the vortex toward the surface is ac- 
companied by the penetration of the antivortex and their 
annihilation, like that in Fig. 10. For higher frequencies, 
this change in the vortex dynamics formally occurs only 
for Bq > B2, so the down step in Q{Bq) does not show 
up in curve 2 in Fig. 12. 



Next we linearize By ~ B^ 
ST,„, where B'^ = drByiTo) 
and B2 take the form: 

S„ = (1 - b)B,„ 

6 = 



— \B'JST,n with respect to 
Then the effective field By 



S2 - (1 - b)B2, 



TT^fiok ln(4(i/7rro) 



(40) 
(41) 



Thus, both B2 and By are reduced by dissipation, which 
only produces constant shifts of local By values but does 
not change their initial distribution. For Nb, taking k — 
10 W/mK, B'y ~ By/Tc, By = 0.15 T, uj/27: = 2 GHz, 
ln(4c?/7rro) ~ 10, as before, we obtain a small value b ~ 
3 X 10^'*, for which the shift of By is negligible. However, 
for NbaSn, with k ^ 10~^ W/mK, we get a much higher 
value b ~ 0.3, indicating that dissipation can significantly 
reduce By and B2, expanding the field region By < Bq < 
Bs of individual vortex penetration. 

The condition of the single-vortex penetration By < 
Bq < B2 implies that the local value of B2 is smaller than 
the uniform superheating field Bg- Multiple vortex pen- 
etration for Bq > B2 causes strong dissipation, further 
decreasing B2 and resulting in avalanche-type dendritic 
vortex penetration™'- 1. Such thermo-magnetic dendritic 
flux avalanches have been observed in both low-Tc and 
high-Tc superconductorsSSiSiHiS^iS^iSl. Notice that the 
superfast vortex penetration through the surface bar- 
rier due to the jumpwise LO vortex instability may per- 
tain to the supersonic vortex velocities observed for den- 
dritic vortex penetration in YBa2Cu307 and YNi2B2C 
films^^S. 

Temperature distribution (|29p results in the long-range 
dc repulsion force iriL) — — s*VT between two oscillat- 
ing vortices spaced by L. Here s*(T) is the vortex trans- 
port entropy responsible for thermomagnetic effects in 
the mixed state®*'^^. In thick films d ^ X, vortices are 
localized at the surface, so to calculate the thermal force 
friL), we put a; = in Eq. (|29p . and obtain 



ML) 



s*Q 



2dk smh{nL/2d) 



(42) 



This long-range force on the scale much greater than 
A results in repulsion of neighboring vortex penetration 
channels, facilitating bending instability and dendritic 
branching of the multi-vortex tracks. 



D. Thermal self-localization of vortex penetration 



IV. EFFECT OF PINNING 



Local temperature increase around oscillating vortex 
reduces both critical fields By = By{Tm) and B2 — 
B2{Tm) as compared to their isothermal values By{TQ) 
and B2{Tq). To evaluate this effect we combine Eq. 
for Q and Eq. PT|) for T,„ and obtain for L0T2 ^ 1: 



5T„ 



2UJ(j)QBy 



TT^/ioA: ln(4(i/7rro) ' 



(39) 



A. Trapping rf vortices at strong fields, Bo ~ By 

At the initial stage of vortex penetration. Bit) ~ By, 
the driving force is much stronger than typical pinning 
forces by materials defects. However, as the vortex moves 
deeper in the sample, the force Fl{u) cx exp(— m/A) sinwi 
decreases exponentially, so pinning becomes more effec- 
tive if the vortex trajectory passes a pin aligned with the 
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FIG. 13: Vortex pinned by a chain of defects near the sur- 
face. The sohd Une shows the equihbrium vortex shape due to 
competition of pinning and the image attraction forces. The 
dashed hues show instantaneous vortex profiles for B{t) — Bo 
and B{t) = —-Bo, between which the vortex hue osciUates. 



place of the vortex entry within a "belt" ~ A wide. In this 
case the vortex can be trapped and stay pinned as the rf 
field changes sign. However, as B{t) reaches —By during 
the negative rf cycle, an antivortex penetrates along the 
same trajectory as the vortex did and annihilates with 
the pinned vortex. 

The power dissipated due to the vortex-antivortex an- 
nihilation can be evaluated from the change of the ther- 
modynamic potential G{u): 



UJ 
TT 



2i?„ 



-tip/A 



(43) 



This quasi-static expression differs from Eq. (|13p by 
the factor 1 — exp(— Up/A), which accounts for the finite 
Meissner currents at the pin, x = Up, and the term due 
to the gain in the pinning energy Uq as the pinned vortex 
annihilates with the antivortex. 



B. Residual surface resistance at Bo ^ By 

Besides dissipation due to vortex penetration and exit, 
there is also dissipation due to rf oscillations of vortices 
already trapped in the samplers. In this section we ad- 
dress rf dissipation and depinning of vortices trapped in 
a superconductor during field cooling, in the limit of very 
low flux density, for which the intervortex interaction is 
negligible. We first consider a single vortex pinned by a 
chain of equidistant defects spaced by £ from each other 
and by d from the surface in the layer of rf fleld penetra- 
tion, as shown in Fig. 13. In this case the equation of 
motion of the vortex line takes the form 



r]QU = eu + 



boBo _ 



Mo A 



sin ojt 



A 



(44) 



where the first term in the r.h.s describes bending stress 
of the vortex line, the prime means differentiation over 



the coordinate y along the surface, and the last term is 
the sum of the elementary pinning forces fp{x,y). For 
£ > A, the dispersive vortex line tension e reduces to the 
vortex self energy e — (j>oBci/ fio per unit length^-. 

As evident from Fig. 13, the magnetic attraction to 
the surface makes the pinned vortex not straight even at 
zero rf field. As a result, there is a minimum trapping 
distance dm, so that only vortices spaced by u > dm 
can be pinned. Vortices spaced by u < dm are unstable 
and annihilate at the surface, since the image attraction 
prevails over pinning. For weak identical pins, dm can be 
evaluated from the force balance equation: 



2t: fiijX^ 



2d„ 



fp 



(45) 



where fp is the maximum elementary pinning force. The 
vortex segments between the pins bow out toward the 
surface, but for vortices in the trapped flux zone x > dm, 
the curvature of u{y) is weak, and the image attraction 
force Fi = {(j)f^/2TT^o\^)Ki{2u/ X) is nearly uniform. In 
this case the equilibrium shape of the vortex segment 
between the pins is determined by the equation euq = Fi 
with uo(±^/2) — d, which yields the parabolic profile: 



uo{y) = d-uo 



1-7^1' (46) 



UOr, 



IBTT/^oeA-^ 



A 



^'j^i(2d„,/A) 
SAluK 



(47) 



Here e — (?!)§(lnK -|- c„)/27r/ioA^ where the constant 
Ci, ~ 0.5 accounts for the vortex core contribution. It 
is convenient to use the effective k w 1.65k defined by 
hiK = h\K + Cy. Eqs. and (|17|) correspond to £ 3> A, 
so the condition that wom <C A is provided by dm > A, 
as follows from Eq. (|45)) . For denser pins, € < A, the 
nonlocal expression for e should be used^l, in which case 
luK in Eq. ((T7)) is to be replaced by ln(f/^). 

To calculate the power Qy dissipated by the pinned 
vortex under the weak (i?o ^ By) rf field, we seek the 
solution of Eq. in the form u{y,t) = uo{y)+Su{y,t). 
Here the Fourier component 5u^{y) = / 5u{y,t)e~'^'^*dt 
of the oscillating vortex displacement 5u{y^ t) satisfies the 
linearized equation 



(48) 



where/^ = (/)oi3oexp(— d/A)/yUoA. InEq. (|48)) we neglect 
the image force oc exp(—2d/ X)6u, which is much smaller 
than the first elastic term in the r.h.s. for weak pinning 
and d> dm defined by Eq. (gSl). The solution of Eq. (|48)) 
with the boundary condition of the fixed vortex ends at 
the pins, Sui^{±£/2), has the form 



cos[(l 



t)nl/^y/£] 



iLurio 

f^T) — 



cos[(l - i)nl^^/2] 



r)o£y2e. 



(49) 
(50) 
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FIG. 14: Frequency dependencies of the residual surface re- 
sistance due to pinned vortices. Curve 1 corresponds to 
Ri for a periodic chain of pins described by Eq. (|54|l . 
and curve 2 corresponds to Ri for the exponential distri- 
bution of pinned vortex segments given by Eq. (|58|) . Here 
7?i(cx)) =0g{e-2''/^)/A2^oa. 



Here we introduced the pinning relaxation time constant 
Tp and the dimensionless frequency flp. The dissipated 
power per unit vortex length, 



where we used the Bardeen-Stephen rjo and took e = 
(j)llnii/2nfioX^ for £ > A or In k ^ \n{i/^) for ^ <C £ < 
A^^. These results show that: 

1. For LUTp <^ 1, the frequency dependence of Ri is sim- 
ilar to the BCS surface resistance, Rbcs oc exp(— A/T) 
at T ^ Tc- However, unlike i?bcs(r), the vortex contri- 
bution Ri remains finite at T — > 0, so trapped vortices 
can contribute to the non-BCS excess surface resistance, 
which has been often observed on many superconduc- 
tors at low temperatures. In this case even a few pinned 
vortices can result in Ri comparable to the exponen- 
tially small Rhcs{T). This scenario was first suggested 
by Rabinowitzl? who modeled pinning by a phenomeno- 
logical spring constant and did not considered the critical 
depinning spacing dm due vortex attraction to the sur- 
face. The account of a more realistic discrete pin struc- 
ture in Fig. 13, and the gradient of the Lorentz force 
changes the frequency dependence of Ri and results in 
new effects discussed below. 

2. Ri increases significantly as the superconductor gets 
dirtier. To evaluate this effect, we make a rough estimate 
(exp(— 2d/A)) ~ exp(— 2dm/A)A/a, which takes into ac- 
count the main contribution to Ri from vortices in the 
critical belt dm < d < dm + X- Taking exp(— 2fim/A) 
from Eq. (I45|) and using the asymptotic expansion 
Ki{z) = (7r/2z)-'-/^ cxp(— 2), we obtain 



2£ 



can be calculated substituting here Eq. (|49 
forward integration then yields 

sinh z + sin z 



(51) 



A straight- 



T^{z) = 1 



:(cosh ; 



(52) 
(53) 



For uJTp ^ 1, Eq. ([52]) gives the frequency-independent 
limit, Q„rn — > inversely proportional to tjq. How- 

ever, for UJTp <C 1, we obtain the quadratic frequency 
dependence, Qi, = f'^rj^i'^uj'^ /2AQe^ , proportional to 7/0- 

The vortex rf power Q^/a — i?Qi?i/2/^Q per unit surface 
area results in the additional surface resistance Rf. 



Ri — ^rw(-\A^), 



(54) 



where a is a mean spacing between pinned vortices and 
(...) means averaging over the vortex positions d in the 
direction perpendicular to the surface. Since these po- 
sitions must satisfy the stability condition d > dm, the 
main contribution to Ri comes from vortices in the crit- 
ical belt dm < d < dm + A where {e^^'^^^) is maximum. 
For LJTp <^ 1, Eq. ([5i)l simplifies to 



Ri 



15a^(j)QPn In^ K 



(56) 



Eq. shows that Ri cx is proportional to the 

trapped flux density, similar to the Bardeen-Stephen flux 
flow resistivity. In the dirty hmit, A ~ Ao(Co/^i)^^^, 
^ ~ (^o^i)^^^j pn oc l/ii, we obtain that the excess 
resistance Ri oc p'^ increases rapidly as the resistivity in- 
creases (here slowly varying logarithmic dependencies of 
Ri on £i are neglected). This behavior of Ri has been 
observed on Nb cavities in which the change of pn at the 
surface was caused by a low-temperature bakin g^^i^^i^^ . 

3. There is a strong dependence of Ri on the size 
£ of pinned vortex segments because shorter segments 
have stiffer spring constants ^ e/£'^ and thus smaller vi- 
bration amplitudes under rf fleld. If pinning centers are 
distributed randomly, Ri {£) should be therefore averaged 
over the distribution of the segment lengths: 



R,{lu) 



R^{ujrio£^l2e)P{£:)d£ 



(57) 



where Ri is given by Eq. ([M]) . and P{£) is a distribution 
function of the segment lengths. For example, P{£) can 
be taken in the form P{£) = £q^ exp(— £/£o) used in the 
Granato-Liicke model of pinned dislocations, where £q is 
the mean segment lengtb^^. In this case, Eqs. ([SH)) . (|5T 
and dill) yield 



R^ 



GOapn In^ k 



(55) 



Rii^p) — 




e--^;'''ruz)dz, (58) 
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where z = f^p = wr^o^o/^^: ^-^d i?oo = 

(j)Q{e~'^'^^^) / X^-qa is the high-frequency hmit of Ri{uj). 
The resulting behavior of Ri{uf) is shown in Fig. 14: 
for low frequencies, the contribution of weakly pinned 
large- length segments makes Ri{io) higher than Ri{ui), 
while for high frequencies, the contribution of strongly 
pinned small segments makes Ri{uj) lower than Ri{Lu). 
The overall behavior of Ri{Lij) resembles the power law 
dependence Ri oc uj^ with /3 ~ 0.5 — 0.7, which has been 
observed on Pb"^^ and Nb^"^ at 0.1-10 GHz. 

To estimate the pinning time constant r^, we first eval- 
uate the mean length of the vortex segment £o ■ This can 
be done from an estimate of the single-vortex pinning 
force balance fp ~ Jc(j)o£o/ ^J'0, which express £o in terms 
of the depinning critical current density Jc. For core pin- 
ning, fp ~ ^B^TT^^/^o, where ( accounts for the change 
in the condensation energy by the pin due to variation of 
6Tc/Tc and the mean-free path^. Hence, 



C0O 



(59) 



Taking A = 40nm and Jc = lO^A/m^ for Nb at T < 
yields ^ lO'^C^A. The pinning relaxation time Tp can 
then be obtained from Eqs. ((50|) and ([59]): 



1287r2^oP«A2C2j2inK 



(60) 



Taking A w ^ w 40nm, Jc = lO^A/m^, and p„ = lO'^l^m 
for Nb, we obtain Tp[s] ^ 10~^C^- From the rf measure- 
ments of pinning relaxation time ~ 10~* s in Nb^l we 
then deduce C 0.1 and £o ^ lO^A ~ 4yLtm. Here Tp is 
rather sensitive to the value of C determined by details 
of the order parameter suppression at the pin. 

It is instructive to compare Q from an oscillat- 
ing vortex with Q due to the rf electric field Ei ~ 
BqcuX exp(— c?/A) induced in the fixed normal vortex core. 
In the latter case the power Qy ~ iv^'^Ef / pn gives the sur- 
face resistance Ri ~ nuP'X'^ eyL-pi—^d/ X) / fjL^K^apm which 
is by the factor k"** In^ k <^ 1 smaller than Ri given by 
Eq. ([55)1 for an oscillating vortex segment. Thus, for 
type-II superconductors with k 3> 1 considered in this 
paper, the inductive contributior>2!£ is negligible. 



Low-field nonlinear surface resistance and rf 
annealing of trapped magnetic flux 



some of pinned vortices at d ~ dm can be pushed out of 
the sample by rf field, which means that dm effectively 
increases as Bq increases. 

To calculate the effect of the rf field on the mean at- 
traction force f{d) between a vortex and the surface, we 
average Eq. (|44p in small vortex vibrations 5u{y, t): 



^lKii2d/X) 
2TrfiQX^ 



xt 



t/2 



dy {sin LotSu) 



(61) 



-e/2 



In addition to the first static term in the r.h.s., / contains 
the rf contribution, in which {...)^ means time averaging 
over the rf period. Using {smujtSu)^^ = Ile{Sui^)/2 and 
Eq. (|49p for the Fourier component Su^^^y), we calculate 
the integral in Eq. ([61]) as follows: 



Re 



72 



2i 



fuRe 
rjQUjilp^'^ 



tan[(l - i)nl^^/2] 



1 + i 



(62) 



Taking the real part of this expression, we reduce Eq. 
(ICTj) to the following equation for dm- 



fp ^ 4'lKi{2dm/X) 
£ 27r/ioA3 



^ n£^Blf^[nl'^] 

2/ioA In k 
sinh z — sin z 
z3(cosh z + cos z) ' 



-2d„/A 



(63) 
(64) 



which defines the critical pinning depth dm as a func- 
tion of the rf amplitude and frequency. For wTp <C 1 
and 2dm > A, we can use ra)(0) — > 1/6, and Ki{x) ~ 
(7r/2x)^/2 exp(— x), in which case Eq. ([S5)) becomes 



Bd, — 



90 
ttX£ 



\hr\^- 

V \ 




Ink 




2r(r!y')_ 





-2d„/A 



1/4 



(65) 

(66) 



Eq. (I65p differs from the static Eq. (I45p by the factor 
{I + Bq/B"^), which becomes essential at rather low fields 
Bo <t: By. Substituting e"^'*'"/^ from Eq. ([65]) into Eq. 
([55]) . we obtain the field dependence of Ri: 



RiiBo) 



R^iO) 



l + B^/Bl 



(67) 



So far we have considered Ri independent of the rf field. 
However, because Ri is mostly determined by pinned vor- 
tices in the critical belt dm < d < dm + A, the excess re- 
sistance Ri may become dependent on weak rf field due 
to an increase of the critical distance dm as Bq increases. 
This effect is evident from Fig. 13, which shows that, be- 
cause both the image attraction force and the Meissner 
rf force increase as the vortex moves closer to the surface 
during the negative rf cycle, the vortex oscillations be- 
come asymmetric and shifted toward the surface. Thus, 



where the zero-field Ri{0) is given by Eq. ([55]) . The 
field B^{ijj) = i?0(O)/\/6r^ increases as uj increases, ap- 
proaching B^{uj) = B^{Q)Vty^ / for ujTp ^ 1. 

The decrease of Ri{Bo) results from the field- induced 
shift of the critical belt dm < d < dm + A away from 
the surface, where the screening of the rf Meissner cur- 
rents reduces vortex dissipation. As a result, the rf field 
irreversibly pumps parallel vortices out of a supercon- 
ductor, resulting in a rf "annealing" of the field-cooled 
trapped magnetic flux. The field dependence of Ri{Bo) 
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FIG. 15: The hysteretic field dependence of Rs{Bo) defined 
by Eq. (IMJ for 7 = 1, = 0.2B,, R,{0) = 0.5Rbcs- The 
descending branch is determined by Eq. (|68p in which Bo 
in the last term stays at the maximum Bq reached on the 
ascending branch. 



is therefore hysteretic: if Bq is first increased to a maxi- 
mum value -Bom and then decreased back to zero, Ri{B) 
on the ascending branch decreases according to Eq. ([S7|) 
and stay equal to i?j(i?om) on the descending branch as 
shown in Fig. 15. This hysteretic behavior enables an 
experimental separation of the vortex contribution to Ri 
from reversible mechanisms due to sound generatior^, 
dielectric losses in the substrate22i^ etc. 

The total surface resistance Kg can be written as a sum 
of the vortex component Ri and the BCS resistance 



Rs = Rbcs 1+7 



R.{0) 



Bl) l + Bl/Bl 



(68) 



where we included the first nonlinear field correction to 
Rbcs due to the Doppler shift of the quasiparticle spec- 
trum, and also pairbreaking and heating effects. In 
the clean limit the Doppler contribution 7^ ^ (A/T)^ 
increases as T decrease o^^i^^i^^'^'^i^^ , while the heating 
component oc Rbcs decreases as T decreases^. Be- 
cause the vortex field B^ is much smaller than Be, the 
total RsiBo) can exhibit a nonmonotonic field depen- 
dence: Rs{Bo) first decreases as Bq increases, reaching 
the minimum Rsm at Bq = Bmin and then increases at 
Bq > Bmin, as shown in Fig. 15. Here Smm can be 
found from Eq. 



-B^iiVi — BcB. 



c^ip] 



jRbcs 



(69) 



The condition that the minimum in Rs{Bq) occurs is 
given by B^^^^^ > 0, that is. 



iRbcs Bl 



(70) 



As illustrated by Fig. 15, the rf field cycling could re- 
duce the surface resistance by irreversibly pumping a 



fraction of trapped fiux out of the sample. The rf flux 
annealing considered in this paper is somewhat analo- 
gous to a directional motion of magnetic flux induced 
by transport ac current, resulting in a dc voltage on a 
superconductor— 



V. NONLINEAR HOTSPOTS IN THE SURFACE 
RESISTANCE 

Localized dissipation due to vortex penetration or os- 
cillation of pinned vortices in thick fllms produces a 
long range temperature distribution, which spreads out 
on the scale ^ 2(i/7r, much greater than A (see Fig. 
8). Even if these temperature variations are weak, 
^r(r) = r(r) - To < To, they can nevertheless pro- 
duce strong variations in the surface resistance Rs of the 
surrounding areas, resulting in nonlinear contributions to 
Rs with very different fleld and frequency dependencies 
than Rbcs{T,LLj). This effect comes from the exponential 
temperature dependence of the BCS surface resistance, 



Rbcs (r) oc — exp 



A 



(5T(r)A' 



rf2 
^0 



(71) 



so that even weak variations 6T(r) < Tq can pro- 
duce strong variations in i?f,cs(r) at low temperatures, 
Tq < i5TA. Indeed, substituting the surface temperature 
distribution ST{0,y) = (Q/tt/c) lncoth(7r2//4£i) from Eq. 
(HH) at a; = into Eq. ^ we obtain: 



Rsiy) 



i?te,(To,c^)coth" (^) , (72) 



a{Bo,To,uj) = Q(Bo,To,c^)A(To)M(To)T2 (73) 

on the scales \y\ > vq greater than the size rg of the heat 
source. Here the exponent a is proportional to the dis- 
sipated power Q, which depends on both Bq and cu. For 
example, Q cx {Bq — B^)""^^ near the onset of the sin- 
gle vortex penetration (see Eq. p4|) . or Q oc BQri^{Qp) 
for a pinned vortex near the surface (see Eq. ([5^ ). In 
turn, the dependencies of a on Bq and to result in a non- 
linear contribution to the global surface resistance from 
sparse "hotspots" of size of the film thickness around 
much smaller heat sources^S,. These hotspots contribu- 
tions can have very different dependencies on Bq and uj as 
compared to the field-independent Rbcs oc ut'^. For a < 1, 
the total excess resistance S'Sis is insensitive to the power 
distribution q{x), and can be obtained by integrating Eq. 
(|72p using a new variable ip — taiih^ (iry /Ad): 

/oc 
[Rs{y) - Rbcs]dy ^ 
-00 



1 - cr 



R. 



bcs 1 



(74) 



where ip{x) = dhiT{x)/dx, and T{x) is the gamma- 
functionS. For a ^ 1, the expression in the brackets 
reduces to 7r^(T/4, thus 



d^s = Trda{BQ,uj,To)Rbc 



cr < 1. 



(75) 
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If Q is due to pinned vortices, the correction (53?^ from 
weak hotspots with cr ^ 1 is quadratic in Bp and pro- 
portional to bj'^ for low frequencies ojTp ^ 1 [see Eq. 
([55)) ]. As CT ^ 1, the function in the brackets in Eq. ([7^ 
diverges logarithmically, indicating that the spatial dis- 
tribution of the power density q{x)^ which cuts off the 
logarithmic divergence in T{x,y) in Eqs. should be 
taken into account. 

In the crossover region a ~ 1, the behavior of 
S^{Bq,uj) is sensitive to the details of the power den- 
sity distribution q{x), but in the limiting case cr 3> 1 the 
main contribution to 53? comes from the hottest region 
near the heat source for which dT{y) is given by Eq. ([30)) . 
Substituting Eq. ([30]) into Eq. we obtain 



■bcs 



dy 



-| a/2 



7r2(y2 



r-o 1 — ) (76) 



where = v^T[{a - l)/2]/r(CT/2) w (27r/CT)i/2 for 
(T 3> 1. The behavior of S^s for a > 1 changes radi- 
cally as compared to u < 1: instead of relatively weak 
power-law dependencies of (53?s on Bq and w for a < 1, 
Eq. ([75)) predicts much stronger exponential field and 
frequency dependencies of Sdig for tr >1, tq <C d. The 
case of strong dissipation cr > 1 can result from vortex 
penetration amplified by grain boundarie o^^i^° , surface 
topography^, local enhancements of Rbcs due to impu- 
rity segregation etc. 

The mean surface resistance Rs averaged over all 
hotspot contributions is given by 



Rs = Rbcs + R^+J2 '55R.(rn)M. 



(77) 



Here the averaged residual resistance Ri results from ei- 
ther pinned vortices or other mechanisms, the last term 
in the r.h.s. is due to the effect of vortex dissipation on 
the BCS resistance, r„ are the coordinates of the sparse 
(thermally noninteracting) hotspots, and A is the sur- 
face area exposed to the rf field. As shown above, 6^s 
can have very different temperature and frequency depen- 
dencies as compared to Rbcs , so the hotspot contribution 
can strongly affect the dependencies of global surface re- 
sistance Rs of Lo and T, particularly at low temperatures, 
where Rbcs is exponentially small. Moreover, the last 2 
terms in the r.h.s. of Eq. (|77p can bring about a strong 
dependence of Rs on the rf amplitude, which can be well 
below the field ^ BcT/Tc of intrinsic nonlinearities of 
the BCS surface resistance due to the Doppler shift of 
quasiparticle energie a^'^'^^i^^'^'^i^^ . 



spacing between hotspots Li is shorter than the thermal 
length Lg = {dk/agY^'^, over which 5T{v) decays away 
from a single hotspot. Here ag = kax / {dax + k) is 
the effective thermal resistance across the film which ac- 
counts for the resistance d/k due to thermal conductivity 
plus the interface thermal resistance X/ax, where aK is 
the Kapitza heat transfer coefRcient. For Li <C Lg, ther- 
mal fields of hotspots overlap and the temperature Tm 
along the surface becomes uniform. In this case the ther- 
mal balance equation takes the form 



{Tm - To)ae 



3. 



Ri + Rbcs{Tm)\, 



(78) 



which determines self-consistently both the rf dissipated 
power and the maximum temperature Tm as functions 
of _Bo and lo. It is convenient to express Bf){Tm) as a 
function of T„, from Eq. ([75)) : 



R, 



2nl{T,n - To)ag 

-Rocxp[{Tm-To)A/TS] 



(79) 



Here we took into account the most essential exponential 
temperature dependence of the BCS surface resistance, 
where Rq = Rbcs{To), and the residual resistance Ri due 
to trapped vortices is assumed temperature-independent 
for T <C Tc- The function BolTm) has a maximum at: 



To 



rp2 
^0 



B^Ri 
2^lag ' 



(80) 



giving the critical overheating Tm ^Ta ^ Tq above which 
a thermal instability develops. From Eqs. ([50)) and ()79p . 
we obtain the equation for the maximum field Bp: 



RpB^eA 
MT-^ag 



exp 



RrBjA 

2filT^ae 



(81) 



The thermal balance Eq. ([75)) has solutions TmiBo) only 
if the rf amplitude is below the breakdown field Bp. For 
-Bo > Bp, the thermal runaway occurs because the heat 
generation grows faster than the heat fiux to the coolant 
as Tm increases. This situation is analogous to combus- 
tion in chemical systemg^^ or thermal quench in semicon- 
ductors, normal metals or superconductors^. 

For Ri Rq, the exponential term in Eq. ([5T)) can be 
neglected and the breakdown field is given by-22. 



Br, 



fJ-o 



f2aeTl 
\RoeA 



1/2 



(82) 



A. Thermal instabilities ignited by hotspots 

In the previous section we considered the rf power as 
a function on the bath temperature Tq. However, Q is 
actually determined by the local temperature Tm, which 
should be calculated self-consistently from the heat bal- 
ance condition. We consider the case, for which the mean 



The temperature dependence of Bp{Tq) can be obtained 
taking i?o(To) = i?„o(A/To) exp(-A/ro), ag = a'TJ, 
whence Bp{To) oc ro^^+'''/^ exp(-A/2ro) reaches a min- 
imum at Tmin = A/(s -|- 3) and increases as Tq de- 
creases below Tmin- However, for lower temperatures, 
RbcsiTo) becomes smaller than Ri in which case the ex- 
ponential term in Eq. (|8ip dominates, and Bp{To) ~ 
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FIG. 16: The temperature dependence of Bp{To) calcu- 
lated from Eq. ^ for Rtc,{To) = i?„o(A/ro) exp(-A/ro), 
Ri/eR^a = 2 X 10-^ ag = a'Tg, s = 3, and B^o = 
2fila'T^+^/R„o- 

fio{2agTQ / RiA)^/'^ cx decreases as Tq decreases. 

The behavior of Bp{To) is shown in Fig. 16. The max- 
imum in Bp[T) at the optimum temperature Tmax sep- 
arates the regimes controlled by the BCS surface resis- 
tance at To > Tmax and by hotspots due to frozen flux 
or other mechanisms of residual resistance at Tq < Tmin- 
For Nb, (A ft! 18K), the optimum temperature in Fig. 16 
corresponds to Tmax — 2K, while for NbaSn (A « 36K), 
Tmax — 4K, if Ri/Rno is the same for both materials. 
The above consideration based on the linear BCS surface 
resistance assumes that the breakdown field Bp is much 
smaller than the field ~ TBc/Tc at which the intrinsic 
nonlinearities in Rs become important. A significant in- 
crease of the isothermal Rs{Bo) due to these nonlineari- 
ties can strongly affect the thermal breakdown^, limiting 
Bp by the thermodynamic critical field Be- 



VI. DISCUSSION 

The results presented above show that the breakdown 
of the Meissner state by strong rf fields involves su- 
personic vortex penetration through the surface barrier 
weakened by defects, the jumpwise LO-type instability 
and high dissipation even for single vortices. Such dissi- 
pation results in thermal retardation effects and hotspots 
igniting the explosive thermal instability due to the expo- 
nential temperature dependence of the surface resistance. 
These effects are precursors for the avalanche vortex pen- 
etration and dendritic thermo-magnetic instabilitie o^°'^^ . 

At the onset of vortex penetration pinning forces are 
much weaker than the driving forces of the rf Meissner 
currents. Yet, as the vortex moves away from the sur- 
face by the distance ~ A, the Lorentz force decreases 
exponentially, so the vortex can be trapped by pinning 
centers. Such vortices trapped in the thin surface layer 
of rf field penetration during breaking through the sur- 



face barrier or field cooling of the sample can result in 
a temperature-independent residual surface resistance. 
However, because pinning centers are distributed ran- 
domly, the rf power dissipated by pinned vortices Q{u) 
varies very strongly because of the exponential sensitiv- 
ity of Q cx exp(— 2m/A) to the vortex position. This effect 
results in hotspots of vortex dissipation, which peaks for 
vortices spaced from the surface by distances close to the 
minimum distance dmi for which pinning forces can pre- 
vent vortex annihilation at the surface at low fields Bq 
By . The field dependence of dm causes rf flux annealing 
in which vortices are irreversibly pushed out from the 
surface layer. This effect results in a nonlinear hysterctic 
dependence of Ri{B) at low fields, Bq Be, which may 
pertain to the puzzling decrease of the surface resistance 
at low fields Bq ^ 3 — 20mT, which has been often ob- 
served on Nb^i^^ and other superconductor o^^'^" . Yet, 
Eq. ([^ describes well the field dependence Ri oc 1/Bq 
observed on Nb cavitiesS^. 

Besides the field and frequency dependencies of Ri[B), 
another manifestation of the vortex pinning mechanism 
is the hysteretic behavior shown in Fig. 14 due to the 
rf annealing of the trapped flux. However, this mecha- 
nism caused by the gradient of the Lorentz force is only 
effective for the vortex segments parallel to the surface 
and does not affect vortex segments perpendicular to the 
surface. The segments of pinned vortices perpendicular 
to the surface generally give a field-independent contri- 
bution to however if these segments belong to the 
vortex semi-loop trapped at the surface, the rf Lorentz 
force gradient acting on the parallel component of the 
loop can eventually push the whole loop toward the sur- 
face where it shrinks and annihilates. In this case the rf 
annealing decreases Rs , eliminating some of the hotspots 
caused by trapped vortices. This is illustrated by Fig. 
17, which shows 3 different type of vortices trapped at 
the surface. Vortex 1 cannot be pushed out by the rf field 
because only a small segment of it '--^ A is exposed to the 
rf Lorentz force, while the rest part is pinned in the bulk. 
Vortex semi-loop 2 can be pushed out by the rf field, as 
discussed above. Vortex 3 has a parallel segment, which 
however cannot annihilate at the surface because it is 
held back by other pinned segments which are beyond 
the surface layer of the rf field penetration. 

Penetration and trapping of even single vortices at 
low temperatures can significantly increase the expo- 
nentially small Rs, which in turn, decreases the ther- 
mal breakdown field Bp. For example, flux trapped 
during field cooling in the Earth magnetic field Bi ~ 
40/iT, corresponds to the mean intervortex spacing a — 
{(po/ Bi^/'^ ~ 7/im. To estimate Ri for such trapped flux, 
we use Eq. (I54p for strong pinning at uJTp < 1, taking 
(exp(— 2(i/A)) ^ A/a and dm ~ A. Hence, 

2 2 

i?, ^ ^i?,(oo), i?,(oo) = ^, (83) 

where Ri{oo) is the high frequency limit of Ri{uj). For 
Nb, taking 5^2 = 400 mT, A = 40 nm, /)„ = lO'^f^m, 
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FIG. 17: Vortices (shown as dashed lines) trapped near the 
surface by pinning centers (black dots). 



and Bi = 40/iT, we obtain Ri{oo) — 2.5/xr2, much higher 
than the typical values Ri ^ lOnfi observed on high- 
purity Nb^3^ poj. NbaSn, with /9„ = 0.2fj,nm, A ~ 90 nm, 
we obtain Ri{oo) ~ 3.9/ifi. As follows from Eq. ([83| . 
pinning reduces Ri by the factor (ti;rj,)^/30, so Ri ~ lOnil 
corresponds to {coTp)^ 0.1, or Tp ~ 0.5 ns for a;/27r = 1 
GHz. By contrast, the pinning time constant Tp ^ 10^^ s 
measured by Pioszyk et ali^ seems to indicate that their 
Nb sample was in the weak pinning limit, (wTp)^ ^ 1, 
for which Ri ~ Ri(oo) is in agreement with the measured 
Ri{0.5Oe) ~ 2fin and the estimate from Eq. ([83]) . 

Introducing dense pinning structures in the surface 
layer of the rf field penetration can therefore impede vor- 
tex oscillations and significantly reduce the part of Ri 
caused by trapped flux, particularly for vortex loops like 
1 and 3 in Fig. 17 which are not affected by the rf flux 
annealing. Because pinning is only effective if toTp < 1, 
decreasing Tp in Eq. (|50p implies reducing the pin spac- 
ing. At the same time, a more effective rf flux annealing 
requires both weak and dense pins (small fp and £ in Eq. 
([55)1 ). Overall, Ri can be reduced by decreasing the re- 
laxation time constant Tp, which can be done not only by 
decreasing the pin spacing £ but also by optimizing the 
mean-free path at the surface. Indeed, Eq. ([50]) shows 
that Tp is larger for higher-K superconductors because of 
the softening of the vortex line tension e, although this 
effect can be offset by a higher normal resistivity. For 
example, the ratio Tpi/Tp2 = K\pn2/ n'ipni oc p„i/p„2 for 
two different materials 1 and 2 (or two different mean free 
paths £i) but the same £, shows that pinning becomes less 
effective for a dirtier surface. Furthermore, comparing 
Nb with Ki = 1, pni = lO^^rim and NbaSn with K2 — 30 
and pni = 0.2/ir2m, we obtain r^Vr^''^^^" - 1/5. Thus, 
reduction of Ri by pinning turns out to be somewhat 
more effective in Nb, although this conclusion can be 
strongly affected by impurities, as discussed above. 

The results of this work show that reducing vortex 
dissipation is an important problem in achieving ulti- 
mate pairbreaking breakdown fields in superconductors. 
In particular, a significant progress has been made in 
increasing Rs and Bp by low-temperature annealing of 
Nb cavities which enables tuning the impurity concentra- 



tion, nanoscale oxide layers and hotspot distribution on 
Nb surface22i2ii2^. Another possibility in raising the ulti- 
mate breakdown fields is to use thin film superconductor- 
insulator-superconducting (SIS) multilayer coating with 
high- Be films of thickness d < X to significantly increase 
Bel and delay the field onset of vortex penetration^i. 
Moreover, the SIS coating may suppress the LO instabil- 
ity by decreasing the vortex flight time through the film 
and providing strong pinning due to magnetic interaction 
of the vortex with the film surfaces. The SIS multilayer 
coating of Nb cavities may enable increasing the ultimate 
breakdown fields Bp above B^'' by taking advantage of 
A15 superconductors2£ or MgB^^M. with Be > B^'' and 
potentially lower Rbcs- 
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APPENDIX A: TEMPERATURE OF A MOVING 
VORTEX 

Eq. ([M]) can be written in the dimensionless form: 

e = W^0-0 + /3{0^)s^{t)f[x ~ u{t), y], (Al) 

where 9 — {T — To)/{Tc — Tq), and the time, coordinates 
and vortex velocity are normalized by the thermal scales 
tg — C/a, Lg = {k/aY^^ and vg — Lg/tg, respectively, 
/3(6l„0 = Vo{0m)v1/{Tc - To)k, and s{t) = v/vg is the 
dimensionless vortex velocity. The Fourier transform of 
Eq. (|Aip results in the following equation for the Fourier 
components, Op{t) = J 9{x,y,t) exp{—ipr)d'^r: 

9p + {1+ p2)0p = /p.g(t)e-'P==«(*) (A2) 

where g{t) = l3[0,n{t)]s'^{t). The solution of Eq. is: 

Spit) = fp / e-(i+P')*'-'P-«(*-*')g(i - t')dt' (A3) 
Jq 

For steady-state vortex oscillations, the rf field was 
turned on at = — oo, and v{t) in Eq. (IA3[) ac- 
counts for all oscillations preceding the time t. For 
/(r)^ = TT-^i^'expi-ryei), and /p = exp(-p2^iV4), 
the inverse Fourier transform of Eq. (|A3|) gives 

, 1 dt'n(t-t') t' ^ 

Jo 4t' + ^2 

Eq. (|A4p was obtained for an infinite sample. For a 
thermally- insulated or ideally cooled surface at a; = 0, 
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Eq. (jA4p can be modified using the method of images: 



e e 



dt'g(t - t') 
W + 12 

[x-^(t-t')]^ [^ + u(t-t')f 

£^+4t' -|- g 5^+4t' 



(A5) 



where the plus and minus signs in the parenthesis corre- 
spond to the Dirichlet [dxd{0, y, t) — 0] and the Neumann 
[0{0, X, t) = 0] boundary conditions at the surface, re- 
spectively. The core size is typically much smaller than 
Lg, so the dynamic equation for the core temperature 
Om{t) can be obtained from the self-consistency condi- 
tion Omit) = 9{u{t),0,t), resulting in Eq. (f25l) . 

Next we consider the steady-state temperature distri- 
bution T(x, y) averaged over high-frequency vortex oscil- 
lations, at 2nu}tg » 1. In this case T{x,y) satisfies the 
static thermal diffusion equation 



kW^T + q{x)S{y) ^ 



where q{x) is the mean power density distribution along 
the vortex trajectory. We solve Eq. (jA6|) for a film of 
thickness d with the boundary conditions: dxT{x, y) = 
on the thermally-insulated surface a: = where vortex 
dissipation is localized, and ^^(a;, y) = Tq at the opposite 
surface, x = d kept at Tq. By symmetry, this geome- 
try has the same T(r) as in a film of thickness 2d with 
isothermal boundary conditions T{±d,y) = Tq, and the 
heat source in the middle at x « 0. In this case Eq. (|A6[) 
can be solved using the Green function 



1 cosh -I- cos 

G r,r' = -^In 

47rfc cosh i^^- cos ^^^^ 

which gives the distribution of ST{r) = T{r) — Tq 



(A7) 



^T(r) = — / ,(.')ln 



cosh II -I- cos 



7t{x-\-x') 

2d 



Try 
2d 



7t{x — x') 

2d 



dx' . 



(A8) 

If the length of dissipation source rg is much smaller 
than the film thickness, dT(r) around the source at 
(x^, y^) <C d"^ reduces to 



16d2 



[y2 + {x — x') 



-^dx'. (A9) 



Next, we take a rectangular approximation q{x) = qo for 
x < ro and q{x) = for x > tq, where rg is a charac- 
teristic size of the dissipation source, so that qoro — Q 
gives the total power Q. In this case the distribution of 
i5r(0, y) along the surface becomes 



(AlO) 



which yields after integration: 
16d^ 



ST(y) = 



In- 



(All) 




FIG. 18: Graphic solution of the equation F{v) — F, where 
F{v) = rjQv/{l + ii^/^o). The branches v± correspond to the 
± signs in Eq. (|34|) . and the dashed lines show the Lorentz 
force F for particular values of u in the stable region (1) and 
the unstable region (2). 



(A6) The maximum 5Tm = 5T{{), 0) is given by: 



^T^.^fmi^ + i 

TTfc \ Trro 



(A12) 



APPENDIX B: MULTIVALUED FRICTION 
FORCE 

Eq. (|^ describes an overdamped vortex driven by 
the force F{u,t) balanced by a nonlinear friction force 
r]v/{l + /vq). The force balance equation has either 
two or no roots, as shown in Fig. 17. The velocity V-{F) 
for the left intersection point vanishes at F = and in- 
creases up to Wo as F increases. The velocity for the 
right intersection point decreases as F increases. 

If < Fm the branch v- (F) describes all smooth parts of 
the vortex trajectory and also provides V-(0) — for the 
initial condition m(0) = 0. For F > Fm, the vortex jumps 
from a; = ui to the point x = U2 where friction is able 
to balance the drive. Here mi_2 are defined by the con- 
dition -F(ui,2, ^1,2) — Fm- The branch V-{F) describes 
all smooth parts of u{t) provided that F{ui^2,t) < 0, so 
that F{t) always decreases below F^ after the jump. 

For high rf frequencies or strong vortex core overheat- 
ing, there are situations, for which F > after the jump. 
For example, if u > A, the term with the K{x) in F can 
be neglected, and F — vodxF + dtF becomes: 



F = {~vo/\ + ujcoiujt)F, 



(Bl) 



which tends to become positive at higher frequencies. 
In this case the vortex cannot jump to the point where 
F{u,t) = Fjn, since the friction force cannot balance F 
if > 0, because F{u,t) keeps increasing above F^. so 
v{t) becomes greater than vq- Thus, for > 0, friction 
can only stop the vortex jump if w > wq, thus we have to 
consider the branch v+{F) as well. We should therefore 
construct the trajectory, for which the vortex first jumps 
at t = ti to the new point x = ui and acquires the veloc- 
ity Vi > Vq. After that v{t) continuously decreases from 
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Vi > vq at t = ti to vo a.tt = to, as described by the first 
order differential equation ii = v^[F{u,t)]. Here uq and 
to are determined by the equations: 

F[uo, to] = F,„, vod^F + dtF = 0, (B2) 

which state that the F should change sign as the vortex 
reaches the maximum friction force at the critical velocity 
u — 0. Indeed, if F changes sign at any point of the 
descending branch v+{F), the vortex cannot not reach 
vq, so v{t) passes through a minimum at some v > vq and 



then starts accelerating continuously. On the other hand 
if the vortex reaches Vq &i F > 0, the jump instability 
occurs. Therefore Eq. (jB2| provides the only way for a 
stable switch from the descending branch of v+{F) for 
ti < t < to, to the ascending branch of v-{F). Once to 
and uo are found from Eq. (jB2p . the coordinate of the 
jump, Ui, can be calculated for the vortex going backward 
in time, taking the initial condition u = uo, and F = F^ 
and then, solving the equation u = [F{to ~ ty\ with the 
"negative time" i = to — ti, from i = to i = to — ti. 
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